Short-range correlations and spin-mode velocities in ultrathin one-dimensional 

conductors 
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In ultrathin wires positioned on high-K dielectric substrates or nearby metallic gates, electrons 
can form strongly correlated one-dimensional fluids already at rather high electron densities. The 
density-density correlation function, charge compressibility, spin susceptibility, and electron specific 
heat of such fluids are calculated analytically. The results are relevant for transport and thermody- 
namics of carbon nanotube field-effect transistors and semiconductor quantum wires. 



There is a long-standing theoretical prediction that 
one-dimensional (ID) electrons do not obey the conven- 
tional Fermi-liquid theory but instead form a Luttinger 
liquid (LL)i whose fundamental degrees of freedom are 
bosonic modes that separately carry electrons' charge 
and spin. Experimental verification of the LL theory has 
proved to be a challenge. There is, however, evidence 
for a LL-specific suppression of tunneling into carbon 
nanotubes (CN)^ and a spin-charge separation in quasi- 
ID organics.'^ Vexing questions remain: (1) Is there a 
way to demonstrate a truly dramatic departure from the 
Fermi-liquid in real ID systems? (Aside from the spe- 
cial case of quantum Hall edge states."*) And if so: (2) 
What theory should one use to describe such a regime 
in the presence of long-range Coulomb interactions? Fi- 
nally, (3) What phenomena can one expect from this 
deeply non-Fermi-liquid state? Below we propose that 
the desired strong coupling regime can be realized if an 
ultrathin wire, e.g., a CN is (i) placed on a high-K di- 
electric substrate or (ii) brought close to a metallic gate. 
In contrast, the traditional prescription for obtaining a 
strongly correlated regime is to lower the electron den- 
sity n thereby increasing the Coulomb coupling constant 
Tg = l/2naB, where as = S^K/m,e^ is the effective Bohr 
radius. Unfortunately, this route quickly runs into the 
problem of localization by disorder, e.g., random charges 
on the substrate. The advantage of our proposal is that 
the Coulomb potential of these charges would be strongly 
screened, whereas interactions among electrons would be 
affected much less, as shown below. 

We begin with the dielectric substrate case. We call 
a wire of radius R ultrathin ii R <^ as- Experimen- 
tally, enormous as /R ~ 10"^ ratios are achievable in de- 
vices that use zigzag CNs^ placed on the SrTiOs sub- 
stratei^ {R ~ Inm, h?/m^e^ ~ 50 nm, k ~ 200). We 
show that in such an ultrathin wire Coulomb correla- 
tions are enhanced by a large parameter C = hi(aB/R) 
and a correlated regime — Coulomb Tonks Gas (CTG) 
— appears in a window >C~* <C <C 1 of low rg- The 
CTG can be defined as the state where on all but ex- 
ponentially large lengthscales, x < — as exp(7r^/2rs), 
electrons behave as a gas of impenetrable but otherwise 
free fermions. At such x the LL theory, being an asymp- 
totic long-wavelength theory, has no predictive power, 
and the alternative method presented below is needed; 



~ 1/im should be achievable in current CN devices 
The CTG owes its name to a certain similarity it enjoys 
with the Tonks-Girardeau gas of ID cold atoms. ^ 

We find that to the leading order in 1/C <C 1, the 
short-range density correlations in the CTG are identi- 
cal to those of a one-component free Fermi gas. The spin 
correlations are the same as in the ID antiferromagnet. 
We show that the CTG possesses a number of properties 
akin to the 1 ID Wigner crystalj^ including a neg- 

ative compressibility, a high spin susceptibility and elec- 
tron specific heat, and also anomalous finite-temperature 
transport and tunnelin^S properties. Thus, although 
the ID Fermi gas, the CTG, and the Wigner crystal are 
not different thermodynamical phases (they all are LLs), 
significant quantitative differences in short-range proper- 
ties of the electron system in these three regimes cause 
sharp crossover changes in the observables, similar to the 
boson caseii 

Let us proceed to the derivation of these results. 
The crucial insight comes from the two-body problem. 
Consider two electrons with a relative momentum q ^ 
\/a = l/2rsaB interacting via a model potential U{x) — 

I n{\x\ + R). We face the following puzzle. If ^ 1, 
we have the inequality / kKv <^ 1, so that the kinetic 
energy greatly exceeds the characteristic Coulomb energy 
{v = hq/m^ is the relative velocity). Naively, one may 
expect that the Coulomb potential should be a small per- 
turbation. But the reflection coefficient computed in the 
first Born approximation is equal to ihi{Rq)/aBq, which 
is large in the range of Ts that corresponds to the CTG. 
To resolve this puzzle one has to separate the effects of 
the tails of the Coulomb potential (large-x) and of its 
sharp increase at the origin (small- a;). Indeed, let us ex- 
amine the Schrodinger equation for the wavefunction (jj 
of the relative motion 



-no 



1 



= 0, 



(1) 



where ^ = qx is the dimensionless separation, r — 1/qaB, 
and a — Rq. We focus on the case r^r^ ^l,a<Cl. 
The general solution of Eq. ^ is given by 

0(0 = A,W,,^i m + la) + B,W-„,^^{-i\i\ - ^a), (2) 

where s = sgn(0 and W^.,y{z) is the Whittaker func- 
tionjii The constants Ag and Bg must be chosen to en- 



2 



sure the continuity of (j> and at ^ = 0. The Whittaker 
function has the asymptotic behavior— 

W,r,i/2{^0 - exp[-i(e/2) + ir\n{0 - ttt], ^ » r. (3) 

Due to the logarithm in Eq. ((JJ, which results from the 
slowly decaying l/|a;|-tail of the interaction potential, the 
scattered states are not exactly plane waves. Thus, if t{q) 
denotes the transmission amplitude, its phase depends 
on the distance x from the origin to the points where 
the scattered states are measured, t{q) oc eyip[—ir\n{xq)]. 
However, in a wide range of x, from the classical turning 
point X ^ r/q ^ o?' /as to an exponentially large distance 
X ~ g~^exp(l/r) ~ aexp(l/rs), this dependence of t{q) 
is very slow. The phase shift accumulated over this entire 
interval of x is small and can be ignored. This is consis- 
tent with Coulomb potential being a small perturbation 
at such X. Note that the point x ~ which corre- 
sponds to the the characteristic interelectron distance a 
in the many-body problem, is safely within the indicated 
range of x. Therefore, for our purposes wc can define 
the transmission coefficient by t{q) = B+/A^ at = 0. 
With this definition, one can show that t{q) is given by 

t{q) = ^exp(-7^r)[(d/da)M/!,,,l/2(-^a)]-^ (4) 

which entails (cf. Ref. ITlh 



iq 



2 \n{Rq) 



iq~c{q)' 



as 



qag > 1. (5) 



Thus, there exists a window of momenta, l/as <C 9 <C 
C/as, where the Coulomb barrier is effectively impene- 
trable (opaque) 1^(9)1 Ij due to the strong backscat- 
tering at an exponentially short approach distance, x ^ 
exp(— gas). We conclude that while the tails of the 
Coulomb potential act as a small perturbation at mo- 
menta q ^ kp = im and distances a; ~ a, which are 
the most relevant for the many-body problem, the strong 
short-range repulsion yields the effective hard-core con- 
straint for the charge dynamics. Therefore, in the first 
approximation the Coulomb potential is equivalent to a 
very thin and high barrier, i.e., to a ^-function of a large 
strength. Equation Q supports this identification be- 
cause up to 0(l/£)-terms, t{q) coincides with the trans- 
mission amplitude for the potential 



U{x) = {h^/ni^)c{kF)5{x). 



(6) 



Note that in the opposite limit ?- <C 1 of a low-energy 
scattering the tails of the Coulomb potential cannot be 
ignored. The point of interest x ^ q~^ resides deeply in- 
side a classically forbidden region of the Coulomb barrier 
where the wavefunction (/) depends exponentially on x. 
This regime will not be important in what follows, but 
for future reference we quote the counterpart of Eq. jSJl, 
\t{q)\ ~ (7r/£) exp(— Tr/gas). Let us now apply the above 
ideas to the analysis of the many-body problem. 

Consider the limiting case R = +0 first. We have 
c = 00 and t{q) = for all q; hence, the ground-state 



wavefunction ^' has a node whenever coordinates Xj of 
any particles coincide, 1 < j < M . Using an argument 
similar to Lieb-Mattis theoremji^ one can show that ^' 
has no other nodes. Thus, it must have the form 

^ = exp[W{xi, . . . , xm)] Y\_ sin y (xQi - xqj) (7) 



Qi>Qj 



X (-1)Q$(si,...,sm) 



(8) 



where Ql through QM are the indices in the ordered list 
of the electron coordinates < xqi < . . . < xqm < L 
(periodic boundary conditions are assumed), (— 1)*^ is the 
parity of the corresponding permutation, Sj is the spin 
of Qjth electron, and $ is the spin part of the wavefunc- 
tion. The factor exp{W) incorporates the effect of weak 
1/x-tails of the interaction. Apart from this, ^ coincides 
with the ground-state of electrons with infinitely strong 
(5-function repulsioufiSi i.e., the gas of impenetrable but 
otherwise free fermions (GIF). Below we focus on corre- 
lation functions for which W is not needed directly. 

Due to the strict impenetrability built into particle 
exchanges are forbidden, and so neither the parity factor 
(— 1)'^ nor spin are dynamical degrees of freedom. (So, 
at i? = -1-0 extra assumptions are needed to fix $). All 
correlation functions are slaved to those of the density 
operator p{x), e.g., the two-point cluster function. 



h{x)=a^{^\p{x)p{0'm-l, 



(9) 



In the GIF (W = 0) h{x) is the same^^ as in the spin- 
polarized Fermi gas with the Fermi momentum kp — tto, 

h{x) = ~sm'^{kpx)/{kpxf, (10) 
h{q) = -a + e{2kp - q)q/{2kpn), (11) 

where x,q^^ ^ L, 9{z) is the step-function, and hence- 
forth the tilde denotes the Fourier transform. 

The tails of the Coulomb potential cause a correction 
Sh{x) to Eq. ((Tni) . From the analysis of the two-body 
problem, we expect that at not too small x, Sh{x) admits 
a diagrammatic expansion in rg. Since the term (— 1)'5<I> 
does not affect the dynamics, this expansion has identi- 
cally the same form as for one-component fermions, so 
that the standard calculation yields, in the leading order 



Sh 



Ts q , '^kp 

' — J. — ' 

TT^ kpn q 



Sh{2kF)+'-^ln'\z\, z 



q < 2kp, 

2kp 



2kp 



^r^k^ 
TT^ n q^ 



q > 2kp. 



(12) 
0, (13) 
(14) 



As one can see, Sh{x) is a small correction to h{x) 
[Eq. ((117)) ] up to an exponentially large distance ~ 
a exp(7r'^/2rs). We show below that these first-order 
results for 5h{x) smoothly match at a; '--^ a;* with 
the asymptotic long-distance behavior of h{x) computed 
from the LL theory, which is supposed to resum the 
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perturbative series to all orders. We conclude that at 
r^a <C a; <C the first-order perturbation theory applies 
and that the charge correlations in the CTG are indeed 
no different from those of free spinless fermions. 

This result sufhes to demonstrate that the compress- 
ibility of the CTG is negative. We define the inverse 
compressibility by = (Pe/ dv? ^ where e(n) is the en- 
ergy density of the system. To make it finite the Hartree 
term (interaction with a neutralizing uniform backround) 
must be subtracted away. This is an important difference 
from the case of short-range interactions. To the leading 
order in r^, e is equal to the kinetic energy of the spin- 
polarized Fermi gas plus the potential energy evaluated 
using Eq. for h{q). This simple calculation gives 



2e^ 

K 



4r., \ TT 



■ 7 



(15) 



which is indeed negative at 1 /£ ^ ^ 1 . Here 7 is the 
Euler constant. "'^^ With our definition of k, its negative 
sign does not imply any instability of the CTG towards, 
e.g., phase separation. The phase separation would cost 
a large Hartree charging energy that would outweigh any 
gain due to negative term»i4 

Let us briefly make a connection with the LL theory. 
The above result for >c enters the LL machinery through 
the charge stiffness parameter 



KM 



(TT/4:)hv 



hn 



U{q) 



(16) 



Unlike the "classical" definitioufi in a Coulomb LLSjiSiii 
Kp depends not only on >c but also on the interaction po- 
tential U{q) ^ — 2(e^/K) In Rq. This is again because the 
total energy cost of the charge build-up is the sum of the 
negative >ir-term (correlation energyi^) and the large pos- 
titive Hartree term (electrostatic energy). Kp shows up 
[through Vc{q) — VF/Kp{q)] in the dispersion uj — Vc{q)q 
of charge mode. It also determines the low-g behavior of 
the density correlation function: 



h{q) + a -Kp{q)q/ (kpn), 



0. 



(17) 



We notice that Eqs. (|12|l and (|17ll match at q ^ l/x* 
and take it as evidence that Eq. ((T7|l applies at g < 1/a;*, 
while Eq. 1)12(1 is valid at q > Thus, the first-order 

perturbation theory is sufficient for computing h(x) at 
X < x^,. The LL theory can also be used to calculate 
the crossover from Eq. (|13|l to a different type of sin- 
gularity in the immediate vicinity of 5 = 2kp. If de- 
sired, one can use this to study in detail how the asymp- 
totic exp[-^r, \n{x/x^)] decaji^JML of the "2/1:^" (2a- 
periodic) oscillations of h{x) is recovered at large x. 

So far, we have discussed the case of an infinitely thin 
wire. If R is finite, ^' has nodes only at the coincident 
positions of the same-spin particles. However, due to the 
opacity of the Coulomb barrier in the CTG, h(x) is per- 
turbed very slightly. A much more important difference 
is that particle exchanges become allowed and the spin 



acquires some dynamics. Thus, it becomes meaningful 
and interesting to determine <&. Since particle exchanges 
are still highly suppressed, only those between nearest 
neighbors are relevant. In the CTG they are determined 
by the two-body transmission amplitude t{q). Recall that 
the orbital part of the scattered wave depends on the to- 
tal spin S = Si -|- S2 of the two colliding particles, 

cj) = e^«/2 ± e-'«/2 + [;(^) _ i](g.|?|/2 ± g^|?|/2)^ (18) 

where the upper (lower) sign is for the singlet (triplet). 
To the order 0{t), exactly the same asymptotic scat- 
tered wave would result from the exchange coupling 
U^g = (S1S2 - l/4)Ua{x) if at aU k such that \t{k)\ < 1 
we have Ua{k) — 2(11^ /im^,)t{k)k + const, as can be read- 
ily verified via Born approximation. Although the form 
of the short-range potential Ua-{x) in the real space is not 
unique, this has no significance to the first order in the 
small parameter t. Therefore, the spin-spin interaction 
is captured by the Hamiltonian 

f dk f 1 \ 

I 7r\ ^kS-k - -jPkP-k Imi(fc)fc, (19) 

m, J 27r \ 4 / 

where the integration is to be done up to the ultraviolet 
cutoff /cmax ^ c{kp) and s^, pk are the harmonics of the 
spin and charge densities at wavevector k. The physical 
idea expressed by Eq. (fll^ is very similar to that behind 
the familiar antiferromagnetic coupling (4i^/J7)(SiSj+i — 
niTii+i/A) in a large-C/ ID Hubbard model. 

Since the spins are slow degrees of freedom, we can 
average H^r over the fast orbital motion to obtain the 
usual ID S = 1/2 Heisenberg model 

= JE,SjS^+i, J = -\ h{k)lmt{k)kdk. (20) 
■' In J 

Therefore, the spin wavefunction $ for the ground 
and excited states arc given by the appropriate Bethe 
ansatze.^" Substituting Eq. into Eq. (|20|l and keeping 
only terms 0{l/C), we obtain 



J 



3 TO* C + \ni 



C 



< < 1. 



(21) 



In the CTG, J <^ Ep = fi^kp/Sni^, as expected. The 
velocity of the spin excitations in ID Heisenberg model 
iE:^S Va = {n/2){Ja/h) and the spin susceptibility per 
unit length is x<7 = (5MB/2)^(2/7r?i?;o-). This implies that 
Xa of the CTG exceeds that of the Fermi gas Xa (where 
Vcr ~ v/2) by the large factor x<y/Xa = v/2va ~ rsC » 1. 

The low-temperature electron specific heat of a CTG 
is determined by the velocities of the charge and the 
spin modeSiiSi and is dominated by the latter, Ce = 
{TT/3)ikBT/hva) - (2/7r2)(K/e2)(r2/£)fcBT. This is 
large compared to Ce in the Fermi gas because of the 
smallness of v^. Obviously, is relevant for many ob- 
servables. To verify that our strategy for finding Va- is 
correct, we made sure that it reproduces the known exact 
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resultsASJ^ for the ID Hubbard model and for electrons 
with the contact interaction 

IV (7r^/3)(nnVm,c), n< c. (22) 

A gated wire. — Properties similar to those of the CTG 
may also be exhibited by a modestly thin ID wire, 
R < as, if, instead of a high-K dielectric, it is posi- 
tioned a small distance D away from a metallic gate. 
In that case one can model the interaction potential by 
U{x) = e^/\x\ - e'^/Vx'^ + AD"^ at |a;| > R. This model 
was studied recently by Hausler et al.^^ who surmised 
that Kp 1/2 and v„ oc / hi{D / R) at low n. In the 
regime a ^ D"^ /as 3> as the validity of these state- 
ments can be examined in a controlled fashion. The in- 
teraction potential is short-range, U{x) oc aX x ^ a, 
and opaque, so the system is in the GIF limit. After 
a straightforward calculation of t{q) one finds that c in 
Eq. lO has to be replaced hy^ 

as f rr(l/8) £(5/8)1 1 

^ = exp I ^ — + I , (23) 

where coefficient A ^ 1 depends on the behavior of U{x) 
ai X R, and T{z) is Euler's gamma functionJii Since 
the interaction potential is now short-range, subtraction 
of the Hartree term is no longer necessary. We redefine 
>f^^ + U{q = 0) and obtain ~ iT^e^aBn. 

Equation l|16() then implies that Kp ~ 1/2; however, the 
spin velocity, which can be found by substituting Eq. I|23() 
into Eq. H22|l . differs from the surmise of Ref. l2Ti 
Experimental manifestations. — The predicted large dif- 
ference of spin and charge mode velocities {vc and v^) 
can be verified by momentum-resolved photoemission?: 



or tunneling.-— Since CN is currently the best candi- 
date for realizing the Coulomb Tonks Gas regime, we 
have to mention here that the bands of a pristine CN 
have a two-fold valley degeneracy. This can be accom- 
modated into the model via effective spin-1/2 operators 
Ti. Each two-body term in Eq. H2()|l should then be 
replaced by J{S^S^+l + l/4)(2T,T,+i + 1/2). The re- 
sultant Ha- remains integrableSi and one finds that now 
Va = (tt^ /6){hn'^ /m^c). Transport is another powerful 
probe of the CTG and Wigner crystal regimes. In the 
temperature window J <^ ksT <^ Ep, with J ^ IK (a 
crude estimate), ballistic conductance G through a CN 
should be unusual. Adopting the theory of Ref. _9 to 
the two- valley case, we find G — e^//i, which is 1/4 of 
the non-interacting result, i.e., the "0.25-anomaly." On 
the other hand, G measured through a CN with resistive 
(tunneling) contacts may display an unusual power-law 
decrease^Pi with T . The negative compressibility would 
reduce the charging energy of finite wires, and so Eq. I|15|) 
for K can be tested by careful Coulomb blockade experi- 
ments on the CN quantum dots. Finally, the strong elec- 
tron correlations can also be observedS^ in imaging mi- 
croscopy: instead of the 4a-periodic Friedel oscillations 
of a Fermi gas, charge density near boundaries of defects 
would oscillate with period a, i.e., four times smaller, in 
the CTG. With further descrease of n, these oscillations 
would become strongly anharmonic as appropriate for a 
pinned Wigner crystal. Detailed quantitative predictions 
for all such measurements will be discussed elsewhere. 
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